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We study the dynamics of density fluctuations in purely diffusive systems 
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function becomes long-ranged. We then analyze this behavior in the frame- 
work of nonequilibrium fluctuating hydrodynamics. 
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I. INTRODUCTION 



Recently, there has been a great interest in the study of long-range spatial correlations 
in non-equilibrium systems. A wide variety of models governed by Langevin-like equations 
have been proposed. A first class of models have been introduced to describe systems, in 
the presence of external noise sources, in which the fields are conserved only in average. 
As an example we could mention models of sand-piles JIJ or diffusion in disordered media 
0- ||, which display long-range order. A second group of models would account for sit- 
uations in which internal noise sources are present and the total values of the fields are 
conserved. Among them one finds turbulence models in fluids @], continuum models for 
interface growing or directed polymers || and some models for driven-diffusive systems [[!]], 
@. In them, although fluctuations have an internal origin, they are not necessarily related 
to dissipation. Whereas in the former long-range behavior is always observed, in the latter 
long-range correlations appear when some intrinsic anisotropy is present. 

Hydrodynamic systems constitute a large class of systems which present also internal 
fluctuations coming from the presence of microscopic degrees of freedom. The existence of a 
fluctuation-dissipation theorem introduces new aspects in the behavior of the correlations. In 
particular, it is shown that the form of the fluctuation-dissipation theorem is determinant in 
order to produce such correlations. In this sense, a preliminary result was found for thermal 
conducting systems @] where the temperature correlation function decays as s -1 , s being 
the distance between the points in consideration (see eq.(|3|)). Such long-range correlations 
have been observed experimentally in liquids under thermal gradients |§. 

Our purpose in this paper is to analyze the origin of long-range correlations in purely 
diffusive systems, described only by one variable: the number density of particles. In these 
systems diffusion is the only mechanism responsible for the maintenance of non-equilibrium 
steady states. Therefore, we are not considering the possibility of applying external forces ( 
as for example electric fields) which could be the origin of long-range correlations even if the 
stationary state is homogeneous 0. The dynamics of fluctuations around nonequilibrium 
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steady states is dictated by nonequilibrium fluctuating hydrodynamics. We can imagine a 
general situation in which the diffusion coefficient may depend on the position. This is what 
happens for example when diffusion takes place in a suspension of particles in the presence 
of hydrodynamic interactions or in an inhomogeneous medium. In the first case the diffusion 
coefficient is not constant due to the local concentration value, taking into account the fact 
that the distance among the particles can hinder their movement, while in the second the 
diffusion coefficient depends on the position through the inherent properties of the medium 
where the particles are diffusing. While in the first system we obtain long-range correlations, 
in the second density fluctuations are delta correlated. 

The paper is organized as follows. In section 2 we introduce the models. We specify the 
steady states and the fluctuation evolution equations. These results enable us, in section 3, 
to determine the static correlation functions. Finally, in section 4, we discuss the origin of 
the long-range terms in the framework of non-equilibrium fluctuating hydrodynamics. 



II. DIFFUSIVE SYSTEMS 

Our purpose in this section is to study the fluctuation dynamics of purely diffusive 
systems, which will be characterized by the number density field, n(r,t). According to 
non-equilibrium thermodynamics |[L0|| , this quantity is governed by the continuity equation 



where J is the number density flux, given by Fick's law 

J = D-Vn (2) 

— * 

with D being the diffusion tensor, which becomes a scalar for isotropic systems. To analyze 
the fluctuation dynamics, we will next study two different models with particular expressions 
for the diffusion coefficient. 
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A. Model A 



As a first model, we will consider diffusion in a system in which the diffusion coeffi- 
cient depends on the local number density, as happens, for example, when hydrodynamic 



interactions among suspended particles are considered jLlJ . 

The differential equation governing the evolution of the number density of particles results 
from the continuity equation and the linear law for the mass flux, introduced previously. 
Since the diffusion coefficient is not constant, this equation reads 

^ = V ■ (D{n)Vn) (3) 



dt 

where we have taken into account the fact that the diffusion coefficient depends on position 
through the local value of the density n. The expressions of both the stationary number 
density profile and diffusion coefficient are obtained from the solution of eq. (|3]) in the 
stationary regime 

V (D s (n s ) ■ Vn s ) = (4) 

We will analyze the stationary solutions of (f|) in the situation in which our system fills the 
region between two parallel plates and a concentration gradient Vo^ = {jilii — n -L/2)l L = 
An/L is kept. In this last expression, n± L / 2 are the values of the number density at the plates. 
If the dependence of D s on n is smooth enough, we can expand both the stationary values 
of the diffusion coefficient and the number density in powers of the parameter ^ (fifr) = e ' 
with no being a characteristic value of the concentration. Up to linear order in e one obtains 

n s (r) = n' + f • V n + - — (r • V ^) 2 (5) 

2 n 

where n' = n — |^(An) 2 , with n = (ul/2 + ^-l/2)/2 being an averaged density value, and 
the expression for the diffusion coefficient in this approximation reads 



A*(r) = D 



n e ^± 
1 — e \o n ■ r 

no no 



(6) 
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Note that for small e, D s remains always positive. Once we have obtained the stationary 
solutions, we proceed to study the dynamics of the fluctuations of the number density, 
5n = n — n s . According to fluctuating hydrodynamics, fluctuations follow an equation 
similar to (0) with stochastic sources J R (f,t). Note that care should be taken because 
the fluctuation dynamics will be affected due to the fact that the diffusion coefficient is a 
fluctuating quantity. After linearization in the expansion parameter ( an analysis of the case 
in which mode-coupling terms are present has been given in [12| for a fluid in stationary 
flow), and introducing the expression (H), one obtains 

<}<) " A)(l - e-)V 2 5n - A)V 2 (~r- V n 5n) - V • J R (7) 



dt n \n 

The random part of the stochastic current satisfies the fluctuation-dissipation theorem 

< JV,*) J R (r,t) >=2k B A s D s {r)5{r-r)5{t-t)l (8) 

OjJ 

where 1 is the unit matrix and the average is taken over a stationary ensemble. Therefore, 



6n describes a gaussian process |13|| . In the above expression, the diffusion coefficient is 

dn I 
9/i I 



precisely (01) and the derivative |^| s must be evaluated in the steady state. For ideal systems 



(§) reduces to 

< J R {f,t) J R {r,t) >=2k B Tn s {f)D s {f)5{f-r)5{t-t')f (9) 

which constitutes a good approximation in the case of dilute systems. 

The formal solution for the number density fluctuations, eq.(^), in Fourier space is given 

by 

5n(k, u) = -iG(k, u)k ■ J R (k, tu) (10) 
where we have defined the Green propagator 

G(k,u) = 7 ^ =; =r- (11) 

-iu + Do ( 1 - e^)k 2 - i^k 2 V n ■ V? 
with Vfc being the gradient with respect to the components of the vector k. 



Note that in the absence of inhomogeneities (e = 0), this propagator reduces to 

G (^ HGo (^ ) = _2__ (12) 

which describes the fluctuation dynamics of diffusion processes when the diffusion coefficient 
is a constant. The propagator G(k,u) can also be expanded in powers of e. To linear order 
one gets 



G(k,u) = G {k,uj) 



1 + e—k 2 G (k, u) + i—G (k, u)k 2 V Q n ■ Vr - 2ie—G (k, u) 2 V n ■ k 
n no no 

(13) 



where we have considered that G(k,u) is an operator. This expression will be used in the 
next section to compute correlation functions. 

B. Model B 

We now consider a different model, in which the diffusion process takes place in an inho- 
mogeneous medium in which the diffusion coefficient may depend on the position according 
to the expression 

D{r) = D Q + V D ■ r (14) 

Here D(r) is the diffusion coefficient, D is a constant and Vo-D is a constant vector that will 
be determined by boundary conditions. Variations of the diffusion coefficient are limited to 
ensure the positive character of such a quantity. The former expression could be interpreted 
as an expansion in the parameter L -^p , with L being a characteristic length of the system, 
in which higher order powers of the parameter have been neglected. 

The one- dimensional microscopic version of this model could be a particle performing 
a random walk in the presence of traps whose distribution depends on position ||14|| . It is 
possible to show that the probability distribution satisfies a Fokker-Planck equation in which 
one may identify the position-dependent diffusion coefficient as a function of the probability 
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distribution of the traps. If this distribution varies slowly with the position, one obtains an 



expression similar to (II) 



As in model A, the number density evolves according to 

°-g = V • D(f)Vn (15) 

We will analyze the stationary solutions of fll5[) subjected to the same boundary conditions 
as in the previous model. To linear order in AD /Do, we get 

n s (r) = h + r ■ Vo^ JT^ if' ^o-D) (16) 



1 AnAD 



where we have defined n = n 

We now proceed to study the dynamics of the number density fluctuations. Note that 
fluctuation dynamics in this model will differ from (0) because now the diffusion coefficient 
does not fluctuate. To linear order in the expansion parameter one obtains 

^ = D V 2 5n + f ■ V DV 2 5n + V D ■ V5n - V • J* (17) 
ot 

where use has been made of eq.(|TJ). Here, J R is again the stochastic flux which satisfies 
the fluctuation-dissipation theorem formulated in eq.([|) since the system is ideal. In it, the 
diffusion coefficient is now given by eq . (|14|) . 

The formal solution for the number density fluctuations in Fourier space also follows 
from eq.(|H]), with the appropriate Green function, which now reads 

G(k,u) = i ^ ^ ^ (18) 

-iuj + D k 2 + ik 2 V D ■ Vjg + iV D ■ k 



This expression differs from (|TT|) because now the diffusion coefficient does not fluctuate. 
In the same way as we have done in model A, the propagator can be expanded in powers of 



AD 
D ■ 
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III. CORRELATION FUNCTIONS 



To compute the number density correlation function, for model A, we will use (|T0D and 
the expression for the fluctuation-dissipation theorem (^), in (k,w) representation. We then 
obtain 



< 8n(k, uj)5n(k , u ) > = —G(k, u)G(k ,u)k ■ k D n' + in'eVon ■ 

5{k + k') (19) 



+ iD eV n ■ V fc - ^7 nV n : V fc V fc 



where, due to the convolution in real space, G(k, u) should be understood as an operator. 



The equal-time correlation function in real space follows after Fourier transforming ([19]) 
according to 

< 5n(f,t),8n(r,t) > = / dk^ \ dk e ik r ■ 

(27r) 8 J-oo J-oo 

{ r duje~ iujt r dJe~ iJt < 5n{k,u)5n{k',u') >) (20) 

Lv/— oo J — oo J 

As we have a finite system, one should have used Fourier series in the spatial direction in 
which the external gradient is applied. This fact means that our expressions will be restricted 
for points r and r such that \r— r | << 2L/n, with L being the size of the system. Otherwise, 



the discrete character of the reciprocal space has to be taken into account fI5| . Performing 
the integrals, we finally arrive at 

< 6n(f,t)8n(^t) >= n s (r)8(r- r J ) + J ( f - J V °^ 2 (21) 

2 d 7r^ _ r(|) no \r — r \ a z 

to linear order in e. In the former expression T(x) is the Euler function, and it has been 
given for dimensions d>2. When d=2 logarithmic divergencies appear. As expected, eq.(^Il) 
reduces to its equilibrium counterpart when the gradient is turned off. 

For the sake of simplicity, the long-range contribution in eq.(^l|) has been obtained 
keeping only the first term in the virial expansion of the derivative dn/dfi, which corresponds 
to the ideal case. Potential interactions between particles would introduce additional terms 
in dn/dfi pT| , but long-range behavior is still present. 
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Following the same line of reasoning, we can also compute the correlation function for 
model B introduced in the previous section. Up linear order in |Vo-D|, we arrive at the 
stationary function 



which does not exhibit long-range correlations. This expression differs from the equilibrium 
result since now the local number density enters into the correlation. Equation (^l|) reduces 
to this expression as well, if a constant diffusion coefficient is considered. 

At this point it is worth comparing the above results with the exact ones obtained for 
a thermal diffusive system J/J. Although both processes are diffusive and consequently, at 
the deterministic level, they are described by the same equations, one should emphasize 
the different behavior between concentration and temperature correlation functions in the 
presence of external gradients. Indeed, while the former does not always exhibit long- 
range behavior, in the latter long-range correlations are always present. The equal-time 
temperature correlation function has been found to be 



ent, ks the Boltzmann's constant, p the density and c v the specific heat at constant volume. 
Note that eqs.(|21~D and (f23"D have the same form, which indicates that the modifications in 
the propagator induced by the inhomogeneities in the diffusion coefficient do not modify the 
long-range behavior of the correlations. 

From expression (|2lJ) it is straightforward to compute the static structure factor, which 
for d=3 is 




(22) 




(23) 



where T s (r) is the stationary linear temperature profile, VqT the external temperature gradi- 



S(k) =n s (k) + 



e |V n| 2 
n k 2 



(24) 



An estimation of the non-local contribution, S n i, of eq.(^), given through the second term 
of its right hand side, relative to its equilibrium value, S eq , follows from 
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Sm(k) = e|V n| 2 
S eq (k) n^k 2 

In fluid systems, the concentration gradient will induce a coupling between concentration 
and velocity fluctuations [16|]. This fact leads to a non-local term, proportional to k~ 4 , 
which is more important than the above mentioned correction. However, if the wave vector 
is parallel to the external concentration gradient this coupling may be avoided [17] . For this 
geometry, S n i given through eq.(E3F), is the only non-local correction to the static structure 



factor. For a colloidal suspension of hard spheres the quantity ( [25]) is proportional to the 



volume fraction [] 1 and to the ratios An/n and (kL) 1 . 



IV. DISCUSSION 

We have shown that purely diffusive systems away from equilibrium may exhibit long- 
range correlations. In order to clarify the origin of such correlations we will first analyze in 
detail the nature of the Langevin sources. 

Our starting point is the expression for the entropy production corresponding to a irre- 
versible vectorial process, which is given as a product of a flux-force pair ||18|| . Denoting the 
flux by J and the force by X, one has 

a = J-X (26) 
from which one formulates the linear law 

J = L-X (27) 

— * 

where L is the matrix of phenomenological coefficients. 

In the framework of fluctuating hydrodynamics, fluctuations are incorporated simply by 
adding stochastic sources to the currents. One has 

J = L-X + J B (28) 

where the stochastic current, denoted in general by J R , introduces a gaussian white noise 
stochastic process with zero mean and whose correlations are given by 
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< Jk(r, t) J^(r , *') >= f<5(r - r )5(t - t') (29) 

The matrix T is related to the matrix of phenomenological coefficients by means of the 
relation T = 2k B L. This last expression, together with eq.(^), constitutes the formulation 
of the fluctuation-dissipation theorem. 

Note that, according to fluctuating hydrodynamics, the former scheme holds when fluc- 
tuations occur around equilibrium or nonequilibrium steady states. In this last case, the 
matrix T may depend on the position since the matrix of phenomenological coefficients may 
depend on local equilibrium quantities to ensure local gaussianity. As a first example, let us 
consider the diffusion of a contaminant. In this case the entropy production reads 

a = J D ■ V| (30) 

from which one derives the linear law 

j D = yi D . (31) 

In these two last expressions, Jd is the diffusion current, \x the chemical potential, T the 
temperature (assumed to be constant) and the corresponding matrix of phenomenological 
coefficients. The linear law ([31]) can be identified with Fick's law, Jd = —D ■ Vn, where 
D = ^Ld^. This identification enables us to formulate the fluctuation-dissipation theorem 

m 



f = 2k B T3^- (32) 
which for non-interacting particles becomes 

f = 2Dn (33) 

One can also apply our former analysis to the case of heat conduction. One arrives at 

f = 2k B T 2 X (34) 
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where A is the thermal conductivity related to the corresponding matrix of phenomenological 
coefficients L q through A = ^L q . 

In the case of heat conduction, for a constant thermal conductivity, the quadratic depen- 
dence of the matrix T on temperature is responsible for the long-range behavior of the static 
correlation function shown in (1231). Concerning mass diffusion, long-range correlations are 
also present only when the corresponding matrix T is quadratic in the thermodynamic field. 
In fact, in model A one has a quadratic dependence of T on n owing to the dependence of the 
diffusion coefficient on the local density field. In model B the diffusion coefficient depends 
on position, but not on the density, and therefore the matrix T is linear in n. 

In order to give a simple justification of the long-range behavior, we may consider the 
following intuitive argument proposed by Ronis ||19|| . If we have a single Langevin source at 
a certain point generating randomly hydrodynamic excitations, it will strongly correlate, at 
equal times, two equidistant points which are closer to this source than the typical hydrody- 
namic decay length. In equilibrium, as we have uncorrelated Langevin sources everywhere, 
all of them being equivalent in intensity and magnitude, the correlations between these two 
points will drop out. Out of equilibrium, however, as the Langevin forces are detuned, 
long-range correlations could in principle appear. 

In the models we have analyzed, while the Langevin forces are detuned long-range cor- 
relations are not always present. Consequently, the previous argument should be somehow 
precised. In model A, and also in thermal diffusive systems, the ratio between the intensity 
of the noise T and the local value of the corresponding field is proportional to the stationary 
field, whereas in model B it does not depend on it. It is our contention that long-range 
behavior originates from the nature of the Langevin sources, which should increase their 
intensity relative to the local thermodynamic field, and not only from their detuned charac- 
ter. Moreover, the form of eq.(|2T|) indicates that the modifications in the propagator due to 
the inhomogeneities in the diffusion coefficient and its fluctuations do not contribute to the 
long-range behavior of the correlations, which appear only through the above-named mech- 
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anism. Finally, note that the correlation functions decay as k~ 2 as we have purely diffusive 
systems and there is no coupling with the fluctuations of other hydrodynamic fields. 

Long-range correlations have also been predicted in other nonequilibrium systems, as for 
example for fluids under temperature gradients |[20|| , for diffusive systems in the presence of 
chemical reactions [T3] or in semiconductors in the presence of an electric field || . In fluid 



systems, correlations between density or temperature and velocity appear out of equilibrium 
owing to the breaking of time reversal symmetry and have been observed experimentally in 
Brillouin scattering What regards Rayleigh scattering, both in simple fluids |§ and in 



binary mixtures JT6[, long-range correlations appear due to the resonant coupling between 
temperature or concentration fluctuations and transverse velocity fluctuations f22|| , which 
leads to the characteristic k~ A behavior of the static structure factor. In the remaining 
examples, long-range behavior is related to the existence of characteristic frequencies or 
lengths. Note that in our case, as in thermal diffusion, detailed balance holds even in the 
presence of external constraints |23[ . 

Contrasting with our analysis, other diffusive systems have been considered in the litera- 
ture, which do not satisfy a fluctuation-dissipation theorem, either because they come from 
an external noise source |TJ or because the internal noise does not have a thermal origin 
@]- |§. In all these examples, the appearance of long-range correlations typically originate 
from some kind of essential anisotropy. However, when fluctuation-dissipation is satisfied, 
anisotropy does not guarantees long-range order. In fact, it is easy to see that model A with 
a constant but anisotropic diffusion matrix does not produce long-range correlations. 

Model A is a real hydrodynamic model in which long-range correlations originate from 
hydrodynamic interactions between particles. A similar model has been proposed as the 
hydrodynamic limit of a lattice gas of interacting particles |[24|| . Though a constant diffusion 
coefficient is considered, the interaction introduces nonlinearities in the number density in 
the fluctuation-dissipation theorem due to the nonideal character of the system. The long- 
range term that is obtained can be deduced from our model A if e is set equal to zero and 
the next order in the virial expansion of dn/dfi is taken into account. 
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